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Abstract Quantum spatial search has been widely studied with most of the study fo¬ 
cusing on quantum walk algorithms. We show that quantum walk algorithms are ex¬ 
tremely sensitive to systematic errors. We present a recursive algorithm which offers 
significant robustness to certain systematic errors. To search N items, our recursive 
algorithm can tolerate errors of size 0{\/^/\nN) which is exponentially better than 
quantum walk algorithms for which tolerable error size is only OiXnN/^/N). Also, 
our algorithm does not need any ancilla qubit. Thus our algorithm is much easier to 
implement experimentally compared to quantum walk algorithms. 

Keywords Quantum spatial search • Robust spatial search • Recursive spatial 
search • Systematic errors 
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1 Introduction 

An important application of quantum algorithms is quantum spatial search (QSS) 
where we search a database of items spatially distributed on the vertices of an under¬ 
lying lattice. The locality constraint must be satisfied which demands that in one time 
step, we can implement only one local operator coupling only neighboring vertices. 
The famous Grover’s algorithm is optimal for general search problems imi2imi4i but 
it becomes inefficient for QSS. However, it can be generalized to get efficient QSS 
algorithms. 

First efficient QSS algorithm was a recursive application of Grover’s algorithm Q 
but subsequent algorithms were based on quantum walk (QW) Ellllll. QW algo¬ 
rithms have been widely investigated both theoretically and experimentally as they 
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offer better time complexity. Unlike their construction and performance, their sensi¬ 
tivity to errors has not been much explored. Errors are inevitable and pose the biggest 
challenge in experimental implementation of algorithms. Completely random errors 
can only be handled using the methods of quantum error correction (QEC) and fault- 
tolerant quantum computation (ETQC). These methods use extra expensive resources 
to add redundancy to the quantum states and gates to overcome small errors. Eor sys¬ 
tematic errors exhibiting specific structures, it is important to investigate the possi¬ 
bility of quantum algorithms which are intrinsically robust to such errors so that we 
don’t need expensive resources used by QEC and ETQC. 

In this letter, first we show that QW algorithms are not robust as they are ex¬ 
tremely sensitive to systematic errors. Then we present a recursive algorithm which 
is intrinsically robust to such errors. We assume systematic errors to be reproducible 
and reversible. Suppose we try to implement an operator 'W but due to errors, what 
we actually implement is iF. Error reproducibility means that every time we try to 
implement we actually implement iF. Error reversibility means that every time 
we try to implement we actually implement . Such errors are not uncommon 
and they arise when there is incorrect calibration of the instrumentation. Eor example, 
they arise due to imperfect pulse calibration and offset effect in NMR systems ii. 

We only consider the case of two-dimensional square lattice as similar ideas can 
be used for other lattice structures and higher dimensions. Our algorithm is based on 
recursion. Unlike the previous recursive algorithm Q, our algorithm does not need 
ancilla qubits to work. Ancilla qubits are a hurdle in the experimental implementation 
and it is preferable to avoid their usage. The original versions of all quantum walk 
(QW) algorithms needed ancilla qubits either to implement coinspace in DTQW or 
spin degrees of freedom in CTQW. But subsequently, it has been shown that QW 
algorithms do not need any ancilla qubit OlOlll llll2lll3lll4lfT5lll6lll7l . a notable work 
being the CTQW algorithm for QSS on a graphene lattice cni which is also experi¬ 
mentally implemented in ifTSll . We extend this advantage to recursive algorithms and 
thus we show that recursive algorithms are easier to implement experimentally. 


2 Robustness of QW algorithms 

Let us first consider the discrete time quantum walk (DTQW) algorithm for QSS 
presented by Ambainis, Kempe and Rivosh (AKR) |0. To analyze its robustness, we 
consider it as a special case of the general quantum search algorithm (GQSA) which 
iterates on |s) to take it close to |f) ifT^ . Here // is the selective phase 

rotation of |f) and Dj can be any unitary operator satisfying Ds\s) = |s). 

Let Ds\£) = be the eigenspectra of with 9g S [—tt, k] and 6^=5 = 0. Let 
0 min < | be the spectral gap and be the two eigenvalues of satisfying 

I'^il < 0miiT Assuming |(s|t)| <C 1 and |A±| <C 0min. the performance of GQSA is 
determined by the moments, 

A = LlWPcott’(0,/2), pG{l,2}. 


(1) 
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Eq. (27) of ifT^ implies that 

where 77 G [ 0 , 7 r/ 2 ] is determined by 

A 


sin277 
Bsin( 0 / 2 ) ’ 


— 


7zB sin 2?] 


= -7ri7,i \P ^ = cot^ +Al, = 1 +A 2 . 

2B|(f|s)| 2 


( 2 ) 

(3) 


The state can be evolved to \ t) using Grover’s algorithm by applying 0(1 / y/P^) 

iterations of Thus the time complexity of GQSA is 

Tgqs = 0 [Pm'^\2q^T[J^] + T[Is] + r[/"])) , (4) 

where T\X\ is the number of time steps needed to implement A. Thus r[A] is 
and7’[AA'] is r[A] + 7’[A']. 

For AKR algorithm, T[L\ is 1 for any local operator L which includes if and Dg. 
But Is is a nonlocal operator coupling each vertex to all other vertices and r[/s] is 
2-JN as shown in IQ. Also, |(f|s)| is \/\/N and 

Take = v^Bsin(0/2) [(;rB/2) +csc 2 t 7] (5) 


In case of no errors, <j) is 7t and Ds is a real orthogonal operator for which A 1 is 0 and 
A 2 is 0(lnA). Thus A is 0, B is 0(-\/lnA) and Take is ©{VNlnN). The systematic 
errors cause disturbance by shifting A from its ideal value of 0, i.e. A ^ 0. For AKR al¬ 
gorithm, cotlrj is A&{^N / InA). If A \nN/\/N then csc2r] ^A©{y^N/\nN) 

7tB/2. Using Eq. (|5]), we get Take = AA0(sin f-). 

A main reason for A 7 ^ 0 is systematic phase errors (SPE) which shift (j) from 
;r to ;r + e. For small errors, e <C 1 and A is e/2 assuming Ai = 0. For e = 2A = 
Q {lnN/\/N), Take is 0 {£)N, much larger than the optimal 0 {Vn) complexity. Typ¬ 
ically, A ^ 1 and In A/ y/N <C 1 hence AKR algorithm is extremely sensitive to SPE. 
This proves analytically the numerical observation of ll20l . It is similar to the sensi¬ 
tivity of search algorithms to SPE due to phase-matching condition 0211122112^ which 
is a source of dominant gate imperfection and poses an intrinsic limitation to the size 
of database A that can be searched m. 

A non-zero Ai also implies A 7 ^ 0. Ideally, Ds is a real orthogonal operator and 
its eigenphase distribution has a symmetry which causes Ai = 0. The existence of Og 
implies the existence of Ogi = —Og such that | (f | 07 )| = |(f |0^')|. Even a slight violation 
of this symmetry can imply |A| = |Ai | ^ inN/y/N and cause AKR algorithm to fail. 
An example was presented in ll^ where Ds is not a perfect real operator. If we know 
Ai then we can choose ^ such that A = cot ^ -f Ai = 0 and AKR algorithm 
works. For the example of ll25l . this was shown in ESI . But (j) can deviate from (j)' 
due to SPE. Thus SPE is a dominant source of error. 

For continuous time quantum walk (CTQW) algorithms, the search is an evolu¬ 
tion under the Hamiltonian H — -f |f) (f |, where is the Laplacian of underlying 
graph of vertices. CTQW algorithms are extremely sensitive to the value of 7 and 
fail for a square lattice if I 7 — 7c| » InA/v^ where jc is the critical value of 7 . This 
is equivalent to the sensitivity of DTQW algorithms to the value of <f>. Thus QW 
algorithms are not robust. Next, we present a robust algorithm for QSS. 
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3 Preliminaries 

Our search space is a x lattice. We assume = 3” for integer n. For 
JCG {1,2,.we define 

Nk- = {0,1,...,3'"-1}, k = n-K. (6) 

We label N lattice vertices by their x and y coordinates where (x,y) G Nn. In the quan¬ 
tum scenario, the vertices (x,y) are encoded by the basis states |x,y) of a V = 3^"- 
dimensional Hilbert space For any K, .^n can be partitioned into 3^’^ = 9”^^^ 
subspaces ^ where (o:,j3) G Nic with each subspace encoding a 3'^ x 3'^ sub¬ 
square. Explicitly, is spanned by the basis states 

\x,y)K^ = \3’"a+x,3'"li+y), (x,y) G Ntc. (7) 

For each subsquare, let 

\sf)=3-’^ E ( 8 ) 

xyGiNic 

be the uniform superposition state (u.s.s.) of all subsquare vertices. 

For any quantum state \x) and angle tJ, let I^{x) be the selective phase rotation 
by angle tt -I- tJ of jj), i.e. 

/"(!) = 1-/#|Z)(X|, U = l-e<-+^'> = l+e‘^. (9) 

If lx) = V|Xo) for a unitary operator V then 

I^{x)=i-UV\xo){Xo\V^=VI^{Xo)vl (10) 

For d- = 0, f'{x) is the selective inversion of \x)- Here tJ is basically an error pa¬ 
rameter. The desired operator is I^{x) corresponding to zero error but due to errors, 
available operator is I^{x)- 

We define the operator as 

n iH^f) = SKiH00K)sl (11) 

a./teNir 

Here5^satisfy5;,|0,0)“^ = |s“^)forall(a,j3)GNK-and/^(00)f) = na,/tG]Njcf^(|0,0)“^). 
We have 

r[/^(i^)]=2r[5J + i, (12) 

where T\X\ is defined earlier and r[/^(00K')] = 1 as /^(OOtf) is a local operator. 

To find we present a method to implement Sk- Suppose the initial state is 

|0,y)”^ for any y. We define the states |/r)a for a G Nk" as 

a :=0 


( 13 ) 
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For a = 0, we choose Lx=o to be a null vector. So |/i)o is the initial state 

We dehne local operators for h G Nff as 

^,Lb\b-l,y)f = \t-l,y)f + ^3’^-l\b,y)f, (14) 

where 3^ = 3'^ — h+1. Using eqs. (I13I14I) and little calculation, we get Li,\ii)i,-i = 
\n)b- We dehne Sk.x as 


= n = Im) 3 --i. ( 15) 

;cGNk- 

As T[Lb] = 1, we have = 3'^ — 1. It is easy to check that |/r) 3 ic_i LxgNk- 

is the normalized state of uniform distribution of amplitudes in .r-direction over the 
vertices of . The roles of x and y directions can be interchanged to design an 
operator 5 k:^ 3 , such that T[ 5 k:, 3 ;] = ^[5k-,v] = 3'^ — 1 and 

S^,y\x,0)f = 3 -’^/^ \x,y)f, (16) 

yGlNlc 

for any value of x. Eqs. (fTsT i and (fThl l imply that the operator Sk = Sk.xSkj satishes 
5^|0,0)“^ = \sf}. Thus 7’[5^] =2r[5^,J =2(3'"- 1) and Eq. (O implies 

r[5^] =2(3'"-!) =4-3'"-3. (17) 


4 Algorithm 

Consider the recursive relation, 

IYk) =UK\t) =I^i\l/K-i)y^{sK)\VfK-i), lv^0) = |f)- ( 18 ) 

Eq. (fTOl) yields I^{yfK) = [1)1]^. Thus 

= Uq = 1n. (19) 

Eq. (fTTl) . with T[I^{t)] = 1 as I^{t) is local, implies 

T[U^] = 3T[U^-i] + T[^\s^)] + T[I^{t)] 

= 37’[t/^_i]+4-3'"-2, T[Uo]=Q. (20) 

Solving above recursive relation, we get 

T[U^] = {4k-1)3’^+1. (21) 


We note that |f) is a unique vertex. Eor any if, let a)3 = a,)3r = T be the index 
of subspace containing the target state |f) and let be its complementary 
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subspace. Let |x,y)^ denote the 9’^ basis states of for (x,y) G Ntf and let \x,y)t'^ 
denote the 9" — 9’^ basis states of . Let |s^) be the u.s.s. of \x^y)\, i.e. 

14) =3-'^ L \x,y)l. (22) 

jtjGNk: 

By definition, C and C . Then Eq. (l22l l implies that the com¬ 
ponent of in each basis state |x,y)^ is Hence = 1/3 or 

l4+i) = (1/3)14)(23) 

We define (Xk = (■s^l Vjc-i) and assume 

(V,c-i|x,y)f = 0 (24) 

for all basis states |x,y)^^. Then leave |</k:-i) unchanged and J^^(sic) = 

/^(s^) for IWk-i)- Thus 

\VK)=I^iWK-i)lH4)\VK-i), (25) 

using Eq. (fTSl l. Then Eq. @ and little calculation imply 

\¥k) = il-fE + fEfs\aK\^)\¥K-i)-a^fs\4)- (26) 

As =^+^1 C Eq. (l24l l implies {y/K-i\x,y)^li = 0. By definition, (i^|x,y)^4 = 
0 and Eq. (l26l l implies that Eq. (l24ll is true for 1 if it is true for any K. It is true for 
tf = 1 as I i/o) = |f) is orthogonal to Hence eq. (l24l i is true for all K. Eq. (|2^ 

implies an immediate consequence Vic-i) = (Xx/'i putting which, along with 

(■S^+I k^) = 1/3. in eq. ([J^, we get 

a^+i = (aj3) [1 -/e-/5+/e/5l«^k] • (27) 

as the recursive relation for a^- 

Eor small errors, we have A = max(|5|, |e|) 1 and fvl— (t?^/2) + it? for 

t? G {5, e}. Retaining only leading second order terms in Eq. dZTl l. we get the follow¬ 
ing relation for (Ok = \cCk\^, 

^ = + + ( 28 ) 

where Wr- = 1 —2cOk- For (0^ < 1/2, > 0 and 1 — (4 cUk:/ 3) is 0(1). Hence cOk > 

(Ok:+i — ©{cOk)- As |i/o) = \t), cti = (■Si I Vo) = 1/3 and ®i = 1/9 < 1/2. Thus, as k 
increases, COk decreases by a constant factor. Suppose, for K = kQ, cOko — 1/lOn. If 
n < kO then CO„> 1 /lOn else Eq. (l28l l implies 

COn > (l/10n)[l-(4«/15)-(2AV3)]''-'"^ 

> 11/150«-AV15, (29) 

retaining only leading order terms and assuming A <C 1 /^/n = 1 /v/hLV. Thus co„ = 
|a„k = f2(l/lnA) as long as errors are small. Using Eq. (fTsT l and Uk = (s^lVic-i). 
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we get |(sn|t/n-i|f)| = = i2(l /y/\nN). This amplitude can be ampli¬ 

fied by using generalized quantum amplitude amplification with arbitrary phases 
5 and e analyzed in ll^ . As long as |e — 5| <C 1 /VlnA^ (which is true as we as¬ 
sume A ^ 1 /VlnA), the amplitude can be amplified to 0(1) by ©{y/\nN) iterations 
of {s„)Un-\I^{t) on Thus the time complexity Q of our algorithm 

is \/\nNT[Un-\\- Putting K = n — I in Eq. (1211 1. we get = ©^/N\nN and 

Q = y/N\r?I^N. 


5 Robustness 

Our algorithm works as long as {5,e} = ©(l/VlnN). This tolerance to SPE is ex¬ 
ponentially better than previous QSS algorithms which fail if e 7 ^ 0(lnN/VN). Due 
to the term in Eq. ([19]) . our algorithm does not just need but 

also their reverse transformations which are available as errors are 

reversible. 

Our algorithm is also robust to local systematic errors in operators S^- If reversible 
errors perturb to then it also perturbs 5]^ to We have assumed S’ to be 
ff-independent for simplicity. As I^{sk) = SkI^ {00k)Sk, eq. (ITST l becomes 

\Wk) = I^{Vk-i)SkIH00k)sUVk-i) 

= lHWK-l)<^^I^i^K)S\WK-l)- (30) 

Defining \Wk-i) =^\Wk-i), above equation becomes {sk)\Wk-i) 

which is similar to Eq. (ITSl) and our analysis holds true provided assumption (i24l) is 
true. By inductive hypothesis, this is true when | y/g) == S\t) is orthogonal to 
Assuming |f) to be the centre vertex of 3 x 3 subsquare , S\t) is orthogonal to 
jS’y for all local errors S. The ideal operator is made up of local operators L\, 
coupling |x,y) either to |x -(- 1 ,y) or |x,y -f 1). It is reasonable to assume the same for 
errors affecting Lj, and hence locality of S is close to real situations. The general case 
of Jf-dependendent S is hard to analyze and beyond the scope of this study. 


6 Discussion 

Our algorithm takes 0(v/iVln^/^A) time steps to search a square lattice. The best 
known performance is 0(v/AlnA) Il27ll28ll^[30l but we need ancilla qubits for this. 
Without ancilla qubits, the best known performance is 0(v/iVln A). The slightly infe¬ 
rior time complexity could be a probable reason why recursive approach to QSS has 
not been as extensively studied as quantum walk based algorithms. 

We have shown that recursive algorithms derive their importance due to signifi¬ 
cant robustness to systematic errors. Our algorithm works as long as errors are of the 
order of 1 /v/lnA and hence error tolerance is exponentially better than that of QW 
algorithms which is InA/v/iV- Also, our algorithm works without any ancilla qubit 
and hence it is easier to implement. These observations indicate that our recursive al¬ 
gorithm for quantum spatial search is more likely to be implementable in near future. 
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Thus it deserves a detailed theoretical and experimental investigation in the same way 
as QW based algorithms are currently being investigated. 
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